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For a set of ten crystalline materials (oxides and semiconductors), we compute the electronic 
band structures using the Tran-Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TB09) functional. 
The band widths and gaps are compared with those from the local-density approximation (LDA) 
functional, many-body perturbation theory (MBPT), and experiments. At the density-functional 
theory (DFT) level, TB09 leads to band gaps in much better agreement with experiments than LDA. 
However, we observe that it globally underestimates, often strongly, the valence (and conduction) 
band widths (more than LDA). MBPT corrections are calculated starting from both LDA and TB09 
eigenenergies and wavefunctions. They lead to a much better agreement with experimental data 
for band widths. The band gaps obtained starting from TB09 are close to those from quasi-particle 
self-consistent GW calculations, at a much reduced cost. Finally, we explore the possibility to tune 
one of the semi-empirical parameters of the TB09 functional in order to obtain simultaneously better 
band gaps and widths. We find that these requirements are conflicting. 

PACS numbers: 71.20.-b, 71.20.Mq, 71.20. Ps, 71.15.-m, 71.15.Mb 



I. INTRODUCTION 

In the last thirty years. Density Functional TheorjD^l 
(DFT) has been successfully applied to the computation 
of structural, mechanical and linear-response properties 
of many different systems.!^ One of the major short- 
comings of DFT is related to the prediction of electronic 
band gaps, known as band gap problem^ Because DFT 
is a ground-state theory, Kohn-Sham eigenvalues cannot 
formally be interpreted as quasi-particle ones. In order 
to obtain reliable band gaps, one should not describe the 
electronic interactions using a mean-field approach such 
as DFT. On the contrary, this is appropriately tackled 
by using Many-Body Perturbation Theory (MBPT), as 
introduced by Hedin and Lundqvist.'^ MBPT is com- 
putationally much more expensive than DFT, and some 
approximations have to be made for real applications. 
Up to now, most of the studies of real systems conducted 
within MBPT have used the so-called one-shot GW ap- 
proximation with a long-standing record of success.'^ 

However, even for the one-shot GW method, the com- 
putational cost is about two orders of magnitude higher 
than standard DFT. Many attempts have therefore been 
made to obtain better band gaps and to improve the 
quality of the excited states without using any computa- 
tionally demanding many-body techniques. The ASCF 
method^ for atoms and molecules, further generalized to 
solids by Chan and Ceder as A-sol,'^ leads to funda- 
mental band gaps in closer agreement with experiment. 
However, the full band structure of the excited states re- 
mains quite approximate. Hybrid functionals can have a 
positive effect on the accuracy of the band gap and the 
position of higher excited state^^^' but are more compu- 



tationally demanding than standard exchange and corre- 
lation (XC) functionals. 

In a recent letter by Tran and Blaha,'^^ a modified ver- 
sion of the Becke- Johnson exchange potentiaP^ combined 
with an LDA correlation part was proposed to improve 
band gap predictions in DFT. This new XC functional 
(referred to as TB09) has already been applied to a large 
variety of system9^2"221 and was able to predict band gaps 
in close agreement with experiment. However, limita- 
tions of this functional have been pointed out by some 
authors. First, there is no energ y functional from which 
the potential is derived jE^I^lEIl which prevents the de- 
scription of any energy-related properties, such as defect 
formation energies or phase stability. Secondly, it is not 
intended to be used for systems without an electronic 
gap and indeed seems to be less accurate for such sys- 
tems.f^Sl F inally , effective masses are reported to be over- 
estimated,^^ with an accompanying narrowing of the 
bands. There are other known drawbacks to the TB09 
functional, namely the fact that it is not size consistent 
and that it is not gauge invariant .^^S! 

In this work, a systematic comparison of electronic 
band structures obtained from the LDA and TB09 func- 
tionals is performed. These are compared with the more 
accurate band structures obtained through perturbative 
GW calculations, and available experimental data. Sev- 
eral oxides and other technologically or theoretically in- 
teresting materials are investigated : silicon, germanium, 
silicon dioxide, tin monoxide and dioxide, zinc oxide, 
calcium sulfide, calcium oxide, magnesium oxide and 
lithium fluoride. We find that TB09 indeed delivers 
band gap in much better agreement with perturbative 
GW calculations and experiment than the simple LDA. 
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By contrast, valence band widths are contracted from 
LDA to TB09 by a narrowing factor varying between 1.0 
and 1.6, and often strongly disagree with both perturba- 
tive GW calculations and experiment.!^ Our perturba- 
tive GW calculations are performed starting from LDA 
eigenvalues and wavefunctions (these will be referred to 
as GoWo@LDA), but also from TB09 eigenvalues and 
wavefunctions (referred to as GoWdQTBOO). The dif- 
ference is non- negligible. In fact, the one-shot GW band 
gaps obtained starting from TB09 are close to those from 
quasi-particle self-consistent GW calculations, at a much 
reduced cost. In order to produce perturbative GW band 
structure, we introduce a new approximate interpolation 
scheme referred to as polyfit. It is based on a piecewise 
polynomial approach, and is especially simple to imple- 
ment. 

In Sec.|lTj the basic theories used throughout this study 
are reviewed. The main ingredients of DFT are first in- 
troduced as well as the approximations to which they 
are intrinsically related. The form of the LDA and TB09 
XC functionals is also specified. Then, implications of 
Hedin's equations and the Many-Body problem are sum- 
marized. In particular, the GW method and the asso- 
ciated technical details are surveyed. Finally, the polyfit 
interpolation scheme is described. In Sec. |III[ the band 
structures, densities of states, band gaps and band widths 
obtained with LDA and TB09 using DFT and GW are 
presented for the ten above-mentioned materials. Then, 
in Sec. |IVt the results and consequences of the use of the 
TB09 XC functional in DFT-based electronic predictions 
are discussed. In particular, we consider the possibility 
to tune the semi-empirical parameter appearing in the 
TB09 in order to achieve better band widths. However, 
this tuning of the band width cannot be achieved with- 
out degrading the description of the band gaps. The final 
conclusions are drawn in Sec. |Vl 



(-2V^ + Vks(r))^P(r) = .P^fS(r), (1) 

where the Kohn-Sham potential VKs(r) relies on the ex- 
ternal potential V'cxt(r) resulting from the electrostatic 
potential of the ions as well as the electrostatic potential 
VH(r) originating from the electronic density n(r) 

VKsir) = Kxt(r) + Vuir) + Kc[n(r)]. (2) 

Therefore, the unique approximation is the XC poten- 
tial, which is usually expressed as a functional derivative 
of the XC energy : V^dnir)] — ^^^^^ . If the exact XC 
functional were known, the method would in principle be 
exact. In practice, only approximations of this functional 
are available for non-trivial systems. There are many pro- 
posed approximations in the literature and the si mpl est 
approach is to use the local density approximatior^J^I for 
the XC energy E^dn] which is expressed as 



E^c[n]= J drn{r)eM, (3) 

with €^c{n) being XC energy density of the homogeneous 
electron gas of density n. In this work, the Perdew- 
Wan^2ll parametrization for the LDA XC energy, based 
on quantum Monte Carlo data from Ceperley and Aldeif^^l 
is used. 

The XC potential proposed by Tran and BlahgP^ is 
composed of a modified version of the Becke-JohnsorP^ 
exchange potential and a LDA correlation part. The ex- 
change part V^^'^^ takes the following form : 



II. METHODS 

All the calculations were done with the ABINIT soft- 
ware package^* using norm-conserving pseudopoten- 
tials. These were prepared using the Atomic Pseudopo- 
tential Engine.'^fil The LiBXC, a library of XC function- 
als,'^^ has been utilized for the evaluation of the XC 
terms. Hereafter, the basic theories, methods and ap- 
proximations that are used throughout this study are 
presented. 



A. Density functional theory 

In DFT, the ground-state of a system is obtained 
by solving a set of self-consistent equations (l][2), as 
proposed by Kohn and Sham.l^l The wavefunctions and 
energies follow from the solution of the one-electron 
Schrodinger equation : 



TB09 



cV^^^\v) + i3c-2) 




(4) 



where V^^{r) is the Becke-Roussel potentiaP^ model- 
ing the Coulomb potential created by the exchange hole 

and ts{r)=^ J2iLi ^"^t ''^i'i is the kinetic energy density. 
The parameter c results from the following equation : 



c — a + 13 



J_ f ^^, |Vn(rO| 



1/2 



(5) 



where Vcdi is the unit cell volume and where a and /3 were 
fitted to experimental gaps in the original paper using a 
least-square procedure. 

In this study, the DFT electronic structures obtained 
using XC potentials from both the LDA and TB09 func- 
tionals are compared. These two sets of KS states are 
then used as initial starting points for perturbative GW 
corrections (GqW^o^LDA and GoH^o@TB09). 
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B. Many-body perturbation theory and the GW 
approximation 

In the Many-body perturbation theory (MBPT), the 
energies and wavefunctions of the system are obtained by 
solving the quasi-particle (QP) equation : 



+ J dr'l](r,r';.. = .?^)Vr(r) = e^^^lr), (6) 

This equation is very similar to Eq. ([T]) except that the 
self-energy E (a non-local, energy-dependent and non- 
hermitian operator) takes the role of the XC potential 
V^c- In a completely general theory, the self-energy is 
obtained by solving a set of five self-consistent integro- 
differential equations w hich were originally introduced 
by Hedin and LundqvistP^ as illustrated schematically 
in Fig. [ija). In practice, the complexity of the problem 
is often reduced by assuming a unitary vertex operator 
r=l leading to the so-called GW approximation as de- 
picted in Fig. [TJb) . The self-energy E is represented as 
the product of the Green's function G and the screened 
Coulomb potential W : 

S(r, r'; ^) - ^ / ^^'^(r, r', u + u')W{r, r', c.') (7) 

It is common practice to avoid the full self-consistency 
by stopping after the first iteration. This has been shown 
to be a good approximation, provided that the initial 
(KS) wavefunctions and eigenvalues are close to the QP 
ones. The resulting method which is known as the one- 
shot GW or GqWq approximation uses a perturbative 
approach for the evaluation of the QP corrections to the 
KS eigenvalues : 



+ Ae. 



GW 



with 



where Z^ ^=1 — (i/;. 



KSI d^ju:) 



\tp^^) is the renormal- 



ization factor and i stands for the spin, A;-point and band 
indexes. 

It is important to stress here that the initial KS start- 
ing point can have a strong influence on the one-shot GW 
electronic structure obtained and in particular, the band 
gaps. For what concerns this study, we highlight that the 
wavefunctions "0^^^ and eigenvalues ef'^ which appear in 
Eq. (|8| strongly depend on the XC functional. 

The Godby and Need^^ plasmon-pole model has been 
used to account for the frequency dependence of the 




FIG. 1. (a) Schematic representation Hedin's equations and 
(b) GW approximation. 



screening in Eq. ([7|. This model has been demonstrated 
to accurately reproduce the states within a large range 
of energies around the Fermi level .122123 Unless explicitly 
mentioned, the extrapolar methocP^l with a compensa- 
tion energy of 2 Ha has been used in order to reduce the 
number of bands included in the computation of the GW 
corrections. 

For eight out of the ten above-mentioned materials 
that we study, we also provide band gaps obtained from 
the Quasi-particle Self-consistent GW (QSGH^) approx- 
imation, introduced by van Schilfgaarde et al!^ Such 
QSGW band gaps do not depend on the starting point 
and tend to slightly overestimate the experimental band 
gap value. 



C. Interpolation of GW eigenvalues : the polyflt 
energy approximation 

The representation of a sufficiently accurate band 
structure (not only the band gap) or density of states 
requires a dense sampling of fc-points in the Brillouin 
zone. It is preferable to avoid such calculations by inter- 
polating the GW eigenvalues on a very dense mesh from 
the results obtained on a coarser fc-point grid. 

One possible technique r elies o n maximally-localized 
Wannier functions (MLWFs).l2alIIlThe production of GI^ 
band structure using Wannier functions is described e.g. 
in Ref. For valence states, this technique easily ap- 
plies. Indeed, the valence wavefunctions are well local- 
ized and the generation of the MLWFs is almost instan- 
taneous. On the other hand, applying this procedure for 
the conduction states is not an easy task as these states 
are usually more delocalized : finding proper Wannier 
functions is far from being automatic. For this reason 
a different and much more straightforward procedure is 
adopted here. 

It is based on the fact that, at least for simple mate- 
rials such as the ones examined here, the GW correction 
to DFT eigenvalues is, to a large extent, only a func- 
tion of the energy. Fluctuations, related to wavevector 
and state dependences, are present, but can be ignored 
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when the computation of the band structure is done for 
plotting purposes. For more complex materials, the GW 
corrections are no longer a simple function of the DFT 
eigenvalues, as observed e.g. in Ref. [331 and the method 
used here may not be sufficient. 

In our polyfit approach, the GW corrections to eigen- 
values are piecewise interpolated using several least- 
squares polynomial fits corresponding to different energy 
ranges. In this way, the dependence of the GW correc- 
tion Aep^ on the band index and A;-point i=(n,k) is 
replaced by a function of the KS eigenvalues only : 



(9) 



First, the electronic bands are split into separate 
groups depending on their energy range, using midpoints 
in forbidden energy gap as limits of the different intervals 
(so-called energy pivots). For each interval, a third-order 
polynomial is employed as a fit for the GW corrections : 



KS 



KS\3 



The polynomial is constrained to go through the end- 
points of each group of bands in order to preserve the 
band widths. Additionally, the spread of the coarse mesh 
eigenvalues with respect to the polynomial is minimized 
by adjusting the polynomial coefficients. Finally, the 
interpolation polynomials are used to correct the DFT 
band structure and DOS on a fine mesh. An example 
of this fitting procedure is given in Fig. [2] for the case 
of a-quartz. The upper panel shows the polynomial fits 
obtained for the GW corrections. 

Such an interpolation procedure is obviously not per- 
fect. The error 9i made for each state i belonging to the 
coarse grid when using this interpolation procedure gives 
an estimate of the overall error of the fit. For each explic- 
itly calculated state on the coarse grid, it is evaluated as 
the difference between the polyfit and the real correction 
Aepw from Eq. M : 



The errors for a-quartz are given in the lower panel 
of Fig. [2] For all the considered valence states and for 
the lowest conduction states (up to 3 eV above the con- 
duction band minimum) , the error is always smaller than 
0.1 eV whereas for higher-energy states displayed in the 
figure, the error goes up to 0.25 eV. For all the systems 
investigated in the present study, the error is less than 
0.2 eV for the valence states and the lowest conduction 
states, and less than 0.5 eV for higher-energy states. 
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FIG. 2. Interpolation of GW eigenvalues with the polyfit en- 
ergy approximation. Example for a-quartz in LDA. Upper 
panel : GoWo corrections for the valence (blue crosses) and 
conduction (red circles) states with respect to DFT eigenval- 
ues. The vertical dashed green lines corresponds to the energy 
pivots chosen to separate the different groups of bands. The 
black lines are the polynomial fits obtained for each group of 
band. Lower panel : Estimation of the error made by the 
polyfit approximation used in the upper panel. 



III. RESULTS 
A. Silicon and germanium 

The experimental silicon lattice parametei^^ 
a=5.43 A was used in this study. Silicon has an 
indirect band gap from F to a point located about 80% 
of the way along the path Y-X. Germanium possesses 
the same zinc-blende crystalline structure (FdSm) with 
an experimental lattice parametei^S a=5.66 A and has 
an indirect band gap from F to L. 

The wavefunctions have been expanded in a plane wave 
basis set with a kinetic energy cutoff of 20 Ha. The k- 
point grids used were 6x6x6 unshifted and 8x8x8 un- 
shifted meshes for silicon (Si) and germanium (Ge) re- 
spectively. For GW calculations, the size of the dielectric 
matrix used was determined by a kinetic energy cutoff of 
10 Ha for Si and 20 Ha for Ge, and 200 bands have been 
used. 

The silicon and germanium band structures and den- 
sities of states (DOS) are presented in Figs. [3] and [4] re- 
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(a) DFT-LDA 

(b) DFT-TB09 



(c) GoWo@LDA 
— (d) GoWo@TB09 
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GW band structures 
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Density of states (a.u.) 



FIG. 3. (Color online) Band structure and DOS of bulk silicon in the diamond structure computed in DFT with (a) the LDA 
XC, (b) the TB09 XC and in one-shot GW using (c) DFT-LDA and (d) DFT-TB09 as a starting point. In the left panel, the 
dashed (black) lines correspond to the DFT-LDA band structure and the full (red/gray) lines represent the DFT-TB09 band 
structure. In the middle panel, the dashed (green/light gray) lines correspond to the GqWqQLDA band structure and the full 
(blue/dark gray) lines represent the GoWo@TB09 band structure. In the right panel, the valence DOS (black line) is compared 
to XPS spectrurrP' in the right panel. 
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(a) DFT-LDA 

(b) DFT-TB09 



(c) GoWqQLDA 
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FIG. 4. (Color online) Band structure and DOS of bulk germanium in the diamond structure computed in DFT with (a) 
the LDA XC, (b) the TB09 XC and in one-shot GW using (c) DFT-LDA and (d) DFT-TB09 as a starting point. The XPS 
spectrum is from Ref. 1461 The color scheme is the same as in Fig. [S] 



spectively. The DOS are compared to experimental X- fair description of the valence DOS. The main features 
ray photoemission spectra"'^ (XPS). It is clear from of the band structure arc indeed correctly reproduced by 
the figures that both LDA and TB09 provide a rather both XC approximations. The more elaborate GqWq- 
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corrected DOS does not give rise to any major change 
in the valence electronic structure. Table |l] gathers the 
direct and indirect band gap values while Table |Il] gath- 
ers the valence band widths as well as first conduction 
band widths, obtained with the different methods. As 
expected, the band gaps obtained with TB09 are fun- 
damentally improved over LDA. GqWo@TB09 opens the 
gap further, overshooting the experimental values. When 
compared to XPS results, our DFT results show that 
both XC functionals yield too small band widths. Fur- 
ther GqWq correction even worsens this narrowing al- 
though the relative error with respect to experiment is 
reasonable in all four cases. 



B. Silicon dioxide 



Silicon dioxide is an in sulato r with a large experimen- 
tal band gap of - 9 eV^2MI ^here DFT in the LDA 
approximation underestimates the gap by ~ 3 eV. It is 
thus a more challenging case than Si and Ge, to test the 
accuracy of the electronic structure obtained with the 
Tran-Blaha functional. Here, the a-quartz polymorph 
(P3221) has been considered with the experimental lat- 
tice parameters and internal coordinates from Wyckofl*^ 
(a=4.91 A, c=5.40 A, usi=0.465, a;o=0.415, ?/o=0.272, 
zo=0.12). An energy cut-off of 40 Ha for plane waves 
and a 4x4x4 /i-point mesh have been used. Finally, the 
dielectric matrix was expanded with an energy cutoff of 
8 Ha and 600 bands were used for the computation of the 
GW corrections. 

Fig. [5] shows the band structures and DOS obtained 
with LDA, TB09, GoWo@LDA and GoW^o@TB09. The 
fundamental gap is indirect, from if to F. The valence 
bands are mainly composed of three groups. The lowest 
bands (not shown in the figure) lie around 20 eV below 
the Fermi level. These bands are mainly composed of 
O 2s-like orbitals. The second lowest group of bands, 
corresponding to bonding O 2p-like orbitals and Si 3s, 3p- 
like orbitals, is located between 5 and 10 eV below the 
Fermi energy. Finally, the upper valence band group is 
composed of non-bonding O 2p-\ike orbitals. The band 
gaps and different band widths of a-quartz obtained with 
the different XC approximations, in DFT and GqWq are 
reported in Table |T] and Table |lT| 

As expected, the TB09 fundamental and direct gaps 
are closer to the experimental values than the LDA ones 
and are similar to recent theoretical gaps^^ using TB09 
in the PAW method. Comparing the DOS with a recent 
XPS experiment,!^ it is clear that the TB09 drastically 
underestimates the valence band widths. Applying GqWq 
corrections rectifies this error, leading to a GoWo@TB09 
valence band structure very close to the GoWo@LDA one. 



C. Zinc oxide 

Zinc oxide (ZnO) is a widely studied material due to 
its application as a transparent conducting oxide. The 
wurtzite structure {PGsmc - BA) of ZnO has been consid- 
ered here with the experimental lattice parameters from 
Kihara and Donnay.^^ It is also a benchmark case for 
band gap predictions as there is no consensus on the the- 
oretical gap obtained from DFT or GW methods, with 
reported GW values ranging from 2.1 eV to 4.2 eV (see 
Ref . l36l and references therein). 

The pseudopotential for the zinc atom includes the 
whole n=3 shell in the valence configuration. A F- 
centered 8x8x5 fc-point mesh has been used and the 
plane- wave energy cut-off used for the wavefunctions and 
dielectric matrix are 150 Ha and 20 Ha respectively. For 
the computation of the GW corrections, 500 bands have 
been included with a compensation energy of 8 Ha for 
higher states. This is equivalent to the use of at least 
3000 bands without the application of the extrapolar 
method. 

The band structures and DOS are given in Fig. [6) The 
DOS is compared to XPS.I^ As already shown in other 
previous studies LDA fails to provide a correct de- 
scription of the Zn 3d levels and gives rise to a strong 
hybridization of these levels with the O and Zn p-states. 
The same authors suggest to add a Hubbard U term in 
order to lower the position of the Zn 3(i-states and ac- 
tually decouple them from the p-states. This results at 
the same time in an increased band gap closer to mea- 
surements. In any case, the same hybridization problem 
is observed when using the TB09 even if a tiny internal 
band gap around -4 eV is obtained. Nevertheless, consid- 
ering all these states as a single group, the valence band 
width is once more clearly narrowed with the TB09 (5.41 
eV) compared to LDA (6.19 eV). The GqWq corrected 
band structures are much closer to each others, at least 
for the highest p states (from -4 eV to the Fermi level). 
In Table [III the present theoretical res ults a re compared 
to the experimental band width (9 eVj^^^ correspond- 
ing to the sum of p states (5.3 eV), d states (2.5 eV) and 
the separating internal band gap (1.2 eV). This exper- 
imental band width is also in good agreement with the 
XPS spectrum reported in Fig. [6j 

The TB09 yields a band gap value of 3.44 eV which is 
again in much better agreement with experimental one 
(3.6 e\S) than the LDA (0.67 eV). The GqWq correction 
to LDA leads to a larger band gap (2.34 eV), still too 
low with respect to the experimental value while starting 
from the TB09 electronic structure leads to a GqWq gap 
of 3.73 eV. 



D. Tin oxides 

Two forms of tin oxides have been considered : stan- 
nous SnO and stannic Sn02 oxide. Stannous oxide crys- 
tallizes in the PA/nmm (BIO) structure while stannic 
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FIG. 5. (Color online) Band structure and density of states of a-quartz computed in DFT with (a) the LDA XC, (b) the TB09 
XC and in one-shot GW using (c) DFT-LDA and (d) DFT-TB09 as a starting point. The XPS spectrum is from Ref.^l The 
color scheme is the same as in Fig. |3] 



oxide adopts the rutile form {PA2/mnm). Their exper- 
imental lattice parameters and internal coordinates are 
taken from Refs. ^BT and US' respectively. For both sys- 
tems, an energy cutoff of 100 Ha was used for the wave- 
functions. The fc-point meshes used were 4x4x3 and 
4x4x6 for SnO and Sn02 respectively. For the GW cor- 
rections, the dielectric matrix was expanded using a cut- 
off energy of 10 Ha and 1000 bands have been used for 
the computation of the GW corrections. 

The band structures and DOS are given in Figs. [7] and 
[8| The highest group of valence bands of both oxides is 
rather well described even if some discrepancies appear 
between LDA and TB09, especially for the band width in 
case of Sn02. For SnO, it is clear that the highly localized 
d states (lying at ~23-24 eV below the Fermi level) are 
better described with the TB09 as they are pushed down 
in energy, closer to the GqWo corrected ones. The O 
s-states (around ~18-20 eV below the Fermi level) are 
better positioned with TB09 than with LDA as confirmed 
by the GgWo results. In contrast with Si02 and ZnO, but 
similarly to Si and Ge, the valence band widths in SnO 



are not narrowed with the TB09 leading to a surprisingly 
good description of the valence states. 

The fundamental band gap of SnO is indirect from F 
to M . TB09 leads to larger fundamental and direct band 
gaps than LDA but both DFT fundamental band gaps 
are still lower than the measured values. The addition of 
a Go Wo leads to indirect band gaps in very close agree- 
ment with the experiment whereas the direct gap at F is 
overestimated by GoWo@TB09. 

For tin dioxide, the flat bands corresponding to d states 
(around -22 cV to -20 cV) are also somewhat better po- 
sitioned in energy with the TB09 than within LDA. The 
group of bands around -16 eV to -18 eV which is com- 
posed of O s states and some Sn s and p states is clearly 
shrunk with the TB09 compared to the LDA and GqWq 
structures. The highest valence bands are also strongly 
contracted with TB09. 

Concerning the band gap, it is clearly underestimated 
within DFT-LDA while the TB09 leads to a value larger 
than the experimental oncI^The GqWq@L\)K band gap 
is still lower than in experiments while the GoWo@TB09 
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FIG. 6. (Color online) Band structure and density of states of ZnO computed in DFT with (a) the LDA XC, (b) the TB09 
XC and in one-shot GW using (c) DFT-LDA and (d) DFT-TB09 as a starting point. The XPS spectrum is from Ref.EH] The 
range of energies from -6 eV to the Fermi level has been magnified with respect to the rest of the XPS spectrum. The color 
scheme is the same as in Fig. [3] 



theoretical gap is closer to the measured value. 



E. Calcium sulfide 

Calcium sulfide has been studied in the rock salt struc- 
ture (FmSm - Bl) using the experimental lattice param- 
eters.f^The energy cut-off used for the wavefunctions was 
45 Ha and the reciprocal space was discretized using 29 
fc-points in the irreducible Brillouin zone. The dielectric 
matrix was expanded using an energy cut-off of 10 Ha 
and 200 bands have been used for the calculation of the 
GW corrections. 

The band structures and DOS are given in Fig.|9] The 
fundamental band gap is indirect from T to X. The TB09 
leads to increased values of the band gaps as shown in 
Table [I] while the LDA underestimates the experiments. 
The GoWo@LDA band gaps are closer to experiment. 
GoWo@TB09 strongly overestimates the direct band gap 
at r while the direct band gap at X and the indirect band 
gap are within 5% of the experimental values. Unfortu- 
nately, no QSGW results were available in literature to 
compare with for this large band gap semiconductor. 

The upper valence band composed of sulfur p-states is 
again clearly narrowed by TB09. In contrast, LDA is in 
much better agreement with the GqWq values. The band 



widths obtained with GqWoQLDA and G'oWo@TB09 are 
in close agreement with each other. 



F. Calcium oxide 

Calcium oxide CaO has been investigated in the rock 
salt structure {Fm3m - Bl). The experimental lattice 
parameter^ have been used. An energy cut-off of 38 
Ha and 29 fc-points in the irreducible Brillouin zone have 
been used for the wavefunctions. The dielectric matrix 
was expanded using an energy cut-off of 10 Ha and 380 
bands have been used for the calculation of the GW cor- 
rections. 

The band structures and DOS obtained are shown in 
Fig. [To] The obtained conduction band structure is some- 
what different with the LDA and the TB09. In LDA, 
there is a local minimum in the conduction band at F 
while this band is pushed upwards with the TB09. In 
both cases, the gap is indirect from F to X. The band 
gaps obtained are gathered in Table |I] There is some con- 
troversy on the exact value and on the indirect character 
of the fundamental gap such that no reliable experimen- 
tal value could be compared with our results. 

Concerning the valence electronic structure, the upper 
valence bands are once more narrowed by TB09 when 
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FIG. 7. (Color online) Band structure and density of states of SnO computed in DFT with (a) the LDA XC, (b) the TB09 
XC and in one-shot GW using (c) DFT-LDA and (d) DFT-TB09 as a starting point. The XPS spectrum is from Ref.Hg] The 
color scheme is the same as in Fig. |3] 



compared to LDA or GoWo-corrected band structures. 



G. Magnesium oxide 

The most stable phase of magnesium oxide (MgO) is 
in the rock-salt structure (FmSm). The experimentaPSl 
lattice parameter a=4.203 Ahas been used. The funda- 
mental gap, which is direct at F, amounts to 7.83 eV.'''* 
The wavefunctions were expanded using 80 Ha as the en- 
ergy cutoff for the plane waves and 44 fc-points in the 
irreducible Brillouin zone. An energy cut-off of 16 Ha 
has been used for the plane wave representation of the 
dielectric matrix while 300 bands have been included for 
the computation of the GW corrections. 

The band structures and DOS are shown in Fig. [TT] 
The band gap is underestimated within DFT-LDA while 
the TB09 leads to a value larger than the experimental 
value. The GoVFoQLDA band gap is closer to exper- 
iments while GoVFo@TB09 pushes the gap even further 
away from the measured value. In fact, the GoT4^@TB09 



gap is very close to the reported QSGW^ result. 

The valence bands are clearly shrunk when using the 
TB09. Comparing the DOS with XPS experimentsP 
shows evidently that the distance between the two peaks 
in the DOS of the upper valence band is much smaller 
than the experimental value when using the TB09 while 
the LDA and the GgWo-corrected band structures are in 
much better agreement. The valence band widths and 
band gap values obtained with the different methods are 
collected in Table Hand Table El 



H. Lithium fluoride 

Lithium fluoride (LiF) has been studied in the rock-salt 
structure [Fniim) using the experimental lattice param- 
eteil^ a=4.028 A. The wavefunctions were expanded us- 
ing 40 Ha as the energy cutoff for the plane waves and 29 
/c-points in the irreducible Brillouin zone. A kinetic en- 
ergy cut-off of 16 Ha has been used for the plane wave rep- 
resentation of the dielectric matrix and 500 bands have 
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FIG. 8. (Color online) Band structure and density of states of Sn02 computed in DFT with (a) the LDA XC, (b) the TB09 
XC and in one-shot GW using (c) DFT-LDA and (d) DFT-TB09 as a starting point. The XPS spectrum is from Ref.Hg] The 
color scheme is the same as in Fig. |3] 



been used for the computation of the GW corrections. 



The band structures and DOS are shown in Fig. 12 



There is a strong narrowing of the valence band when 
using the TB09. Comparing the DOS with XPS experi- 
mcnt^SS] shows obviously that the valence band is much 
smaller than the experimental value when using the TBG9 
while the LDA and the GoWo-corrected band structures 
are in better agreement. The two-peak structure in the 
XPS spectrum is also much better rendered with the LDA 
DOS whereas the TB09 is far from the experiment. Both 
GqWo-DOS fairly reproduce this peak structure. 

The valence band width and band gap values obtained 
with the different methods are gathered in Table |T] and 
Table As for the other materials, the band gap is 
underestimated within DFT-LDA while the TB09 cor- 
rects this, yielding a value of 14.31 eV in much closer 
agreement with the experimental valu^^of 14.2 eV. The 
GqWo correction opens the LDA band gap closer to ex- 
periment while GoH^o@TBG9 opens the gap even further, 
overshooting the experimental value by almost leV. It is 
again very close to the QSGW band gap. 



IV. DISCUSSION 



A. Band gaps 



Fig. 13 shows a comparison of all the theoretical band 
gaps calculated within this study with the experimental 
band gaps. It is clear that TB09 predicts band gaps in 
much better agreement with experiment than LDA. This 
is not surprising as this XC functional was especially de- 
signed for that purpose. Except for CaS and CaO, for 
which the values of the band gaps are not well repro- 
duced, the error is on the order of 5 to 10%. 

The band gaps obtained with GoWo@LDA are in bet- 
ter agreement with experiments than their DFT counter- 
parts. On the other hand, GoWo@TB09 leads to an over- 
estimation of the band gaps for all the compounds of the 
present study. Table |T] also compares our GqWq results to 
published values using Quasi-particle Self-consistent GW 
(QSGW). Both QSGW^ and GoW^o@TB09 tend to over- 
estimate the experimental band gaps. This is in line with 
the usual trend of QSGW to overestimate band gaps. 
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FIG. 9. (Color online) Band structure and density of states of CaS computed in DFT with (a) the LDA XC, (b) the TB09 XC 
and in one-shot GW using (c) DFT-LDA and (d) DFT-TB09 as a starting point. The color scheme is the same as in Fig. [s] 




FIG. 10. (Color online) Band structure and density of states of CaO computed in DFT with (a) the LDA XC, (b) the TB09 
XC and in one-shot GW using (c) DFT-LDA and (d) DFT-TB09 as a starting point. The color scheme is the same as in Fig.js] 



12 



> 



16 
12 
8 
4 




a -4 



-12 
-16 
-20 







n\ 




- 










\ \ 
















(a) DFT-LDA 








— (b) DFT-TB09 
































■ (c)G„Wo@ 


LDA 




- (d)GoWo@ 


TB09 



_1_ 




X W T L W 
DFT band structures 



rr X W T L W 
GW band structures 



r (a) (b) (c) (d) 
Density of states (a.u.) 



FIG. 11. (Color online) Band structure and density of states of MgO computed in DFT with (a) the LDA XC, (b) the TB09 
XC and in one-shot GW using (c) DFT-LDA and (d) DFT-TB09 as a starting point. The XPS spectrum is from Ref.[77| The 
color scheme is the same as in Fig. |3] 



In this respect, GoM^o@TB09 results thus seem closer to 
QSGW ones. To confirm this, one would have to explic- 
itly compare the quasi-particle and TB09 wavefunctions 
but this is beyond the scope of this work. 



B. Band widths 

Conversely, the occupied band widths of most com- 
pounds are contracted when using the TB09 instead of 
the LDA (except for materials with small band gaps, 
namely Si, Ge and SnO). Comparisons with XPS spectra 
and experimental valence band widths clearly indicate 
that this characteristic of the band structure is not de- 



scribed correctly by the TB09. See for instance Fig. 14 
which compares all theoretical band widths to experi- 
mental band widths. As a matter of fact, the TB09 leads 
to a stronger localization than LDA, the latter thus being 
more suitable than TB09 to predict band widths at the 
DFT level. Indeed, band widths obtained with LDA are 
in closer agreement with both experimental and more 
elaborate GqWq results. Finally, the valence electronic 
structures obtained with GoW^oQLDA and Goiyo@TB09 
usually present similar electronic dispersions suggesting 
that the band narrowing obtained with TB09 tends to be 
corrected when using GqWq. 

Defining the narrowing factor as the ratio of the upper 
valence band width in LDA to the upper valence band 
width in TB09, the experimental fundamental gap of the 
materials shows an overall linear correlation to the nar- 
rowing factor as shown in the upper panel of Fig. |15| It 
is indeed clear that the narrowing of these bands is more 
pronounced for materials with wider band gaps. This 



gives some hint on how the TB09 leads to larger gaps. 
Indeed, exact exchange methods such as Hartree-Fock 
lead to an overestimation of the band gaps. Similarly 
the TB09 opens the gaps by mixing in more exchange in 
the system thanks to the adjustment of the c parameter. 
As a consequence of a more attractive exchange poten- 
tial, the electrons are more localized and their electronic 
bands are narrowed. 



Concerning the unoccupied bands, one possible mea- 
sure of their dispersion behavior is the width of the first 
conduction band (ICBW). Such a quantity (the energy 
difference between the minimum and the maximum over 
the Brillouin zone of the lowest conduction band) might 
in principle be obtained from angle resolved inverse pho- 
toemission spectroscopy. However, in this work, we use 
this quantity to analyze theoretical results. Similarly to 
the valence bands, the values for the ICBW reported in 
Table [ll] suggest a narrowing of the unoccupied bands 
within TB09 though not as systematic as for the valence 
bands, as shown in the lower panel of Fig. |15[ 

There has been a lot of interest in this XC functional 
and the TB09 has so far been applied to many systems. 
These studies have been conducted without checking that 
the ground state electronic structure results are coher- 
ent with more accurate MBPT approaches as done in 
this study. Because (a) the band widths are narrowed 
and (b) there is no energy functional from which the 
potential derives, this calls into question the reliability 
of other properties derived from the TB09-DFT ground- 
state. For instance, defect energy levels within the band 
gap or in the valence states might be misplaced. 




FIG. 12. (Color online) Band structure and density of states of LiF computed in DFT with (a) the LDA XC, (b) the TB09 
XC and in one-shot GW using (c) DFT-LDA and (d) DFT-TB09 as a starting point. The XPS spectrum is from Ref.[Hn| The 
color scheme is the same as in Fig. |3] 



C. Tuning the semi-empirical parameter of the 
TB09 functional 

Recently the authors of the TB09 proposed an im- 
provemenlP2i of the original XC by redefining the param- 
eters a and (3 entering in the evaluation of the c parame- 
ter of the functional using a larger set of materials. They 
claim that the TB09 and its revised version yield very 
accurate electronic band structures and gaps. According 
to the present study, the DFT valence band structures 
obtained with the TB09 are not guaranteed to reproduce 
the GqWq band structures nor the experiment. In par- 
ticular, for large band gap materials, the valence band 
widths are strongly contracted. Some studies are con- 
cerned with the determination of the best c parameter 
which can always be tuned such that the experimental 
band gap is obtained. As already noted in the original 
TB09 publication, the band gap increases as the c pa- 
rameter is increased. 

We have examined the possibility to tune the c pa- 
rameter for the materials in which the band width was 
clearly underestimated. Fig. [T6| illustrates this procedure 
for Q;-Si02, ZnO, Sn02 and MgO. Similar trends are also 
observed for the other compounds. The modification of 
the upper valence band width and the conduction band 
is shown for various values of the c parameter ranging 
from 1.0 to 1.7. It is clear that the band widths are 



decreased when the c parameter is increased. This con- 
firms our conjecture that the TB09 leads to band gaps 
larger than LDA by introducing a larger contribution of 
local exchange which at the same time contracts the elec- 
tronic bands. So, no value of this c parameter giving 
reasonable band gaps would correct the TB09's failure 
to describe the band widths correctly. The band widths 
obtained with LDA (black horizontal dashed lines) are 
always closer to the more reliable GqWo results (green 
horizontal continuous lines) than those obtained with the 
TB09. 



V. CONCLUSION 

A systematic investigation of electronic band struc- 
tures obtained either with the LDA or the TB09 XC po- 
tentials has been performed for several oxides and other 
semiconducting and transparent materials. This recently 
proposed functional has attracted considerable attention 
as it is claimed to yield accurate electronic structure and 
band gaps. While the values of the band gaps obtained 
with the TB09 are indeed much improved over LDA, we 
have shown that the overall electronic structure can be 
unreliable especially for large band gap materials. 

Comparisons with GqWq results have been carried out 
for all the materials revealing that the TB09 can strongly 
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TABLE I. Fundamental (Eg) and direct (Eg,d) band gaps of 
all the materials in this study. 
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underestimate the band widths. This contraction of the 
bands is found to be more pronounced for materials with 
a larger band gap. This follows from the fact that the 
TB09 opens the gaps by virtually adding a larger part 
of local exchange to the system through the adjustment 
of the c parameter. This stronger exchange part results 
in a stronger localization of electron hence narrowing the 
band widths. 

As a consequence, the ground-state electronic band 
structure is poorly described for wide band gap mate- 
rials. This is expected to be detrimental for any further 
investigations such as defect studies, analysis of surfaces 
or interfaces, etc. It is shown that any modification of 
the c parameter aiming to reproduce exactly the exper- 
imental band gap would not lead to any improvement. 
Indeed, the band widths of the materials are roughly 
proportional to the inverse of this parameter and will 
invariably be compressed. 
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FIG. 13. (Color online) Comparison of the theoretical and 
experimental band gaps (in eV) for all the different materi- 
als considered in this study. The different XC approxima- 
tions for both DFT and GW calculations are given as black 
dots (DFT-LDA), red squares (DFT-TB09), green diamonds 
(Goiyo@LDA), blue triangles (Go Wo@TB09), and yellow dots 
(QSGiy). Corresponding band gap values are also reported 
in Table |l] Inset: Table of the mean error (ME), the mean 
absolute error (MAE), the mean relative error (MRE, in %), 
and the mean absolute relative error (MARE, in %) for the 
various XC functionals. 



Although very appealing for its efficiency and simplic- 
ity and for its performance in reproducing experimental 
band gaps of pure bulk materials, the TB09 must be used 
with great caution if any other properties than the band 
gap itself are under investigation. The perspectives of 
this kind of approaches are yet promising and there are 
possible avenues for the design of better XC functionals. 
In order to better reflect the actual many-body effects, a 
local XC functional should take into account some part 
of semi-locality such as what is done with the TB09 func- 
tional. 
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TABLE II. Valence band widths (VBW) and first conduction 
band widths (ICBW) of all the materials in this study. See 
the text for the definition of the ICBW. 
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FIG. 14. (Color online) Comparison of the theoretical and 
experimental band widths (in eV) for all the different mate- 
rials considered in this study. The different XC approxima- 
tions for both DFT and GW calculations are given as black 
dots (DFT-LDA), red squares (DFT-TB09), green diamonds 
(GoM^o@LDA), and blue triangles (GoWoQTBOg). Corre- 
sponding band width values are also reported in Table |ll] 
Inset: Table of the mean error (ME), the mean absolute er- 
ror (MAE), the mean relative error (MRE, in %), and the 
mean absolute relative error (MARE, in %) for the various 
XC functionals. 
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FIG. 15. (Color online) Graph of the narrowtng factor of the 
upper valence band (VBWlda/VBWtboq) and of the first 
conduction band (1CBWlda/1CBWtbo9) as a function of 
the fundamental band gap Eg of the materials considered in 
the present study. Materials with larger band gaps exhibit a 
stronger narrowing of the valence bands. 
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